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Two  phenomena  which  are  emitted  from  a  nuclear  weapon,  whose  target 
vulnerability  ranges  often  overlap,  are  thermal  radiation  and  blast  waves. 

In  most  situations  of  interest  to  the  Army,  targets  are  exposed  to  a  high 
percentage  of  the  total  available  thermal  pulse  prior  to  blast  wave  arrival. 
This  sequence  of  loads  is  detrimental  to  the  survival  of  light  weight 
structures  xn  which  significant  in-depth  heating  occurs  prior  to  the  blast 
wave  envelopment.  There  has  been  evidence  that  this  thermal  preheating 
causes  enhanced  structural  damage  [1,2]._~ 

A) 

EXPERIMENT  DESCRIPTION 


In  order  to  demonstrate  these  types  of  synergistic  responses  and  to 
provide  some  quantitative  data,  three  experiments  have  been  conducted  on 
hollow  cylinders  of  aluminum  6061-T6  in  the  BRL  Thermal/”last  Simulator. 

One  experiment  was  conducted  with  thermal-only  exposure,  one  with  blast- 
only  loading,  and  one  with  combined  thermal/blast  loading.  A  diagram  of 
the  Simulator  is  shown  in  Figure  1.  The  target  cylinder  is  clamped  at  each 
end  and  placed  with  its  axis  in  a  vertical  position  as  shown.  The  thermal 
source  is  situated  downstream  from  the  target  to  minimize  the  effects  of 
the  shock  passing  through  heated  air  and  combustion  products.  In  the  com¬ 
bined  thermal/blast  test,  the  target  is  rotated  through  180°  so  that  the 
heated  side  of  the  target  cylinder  is  exposed  to  the  on-coming  shock  wave. 


The  thermal  heating  of  the  aluminum  cylinder  is  accomplished  by  burn¬ 
ing  of  a  mixture  of  aluminum  powder  and  oxygen.  Two  thin  plastic  membranes 
serve  to  contain  the  oxygen  gas  before  ignition.  The  aluminum  powder  is 
sprayed  through  nozzles  into  the  oxygen  atmosphere  and  ignited  starting 
with  the  nozzles  farthest  from  the  target.  The  aluminum  oxide  particles 
resulting  from  the  burn  irradiated  the  near  side  of  the  cylinder,  raising 
its  temperature  by  approximately  200°C. 


The  shock  loading  on  the  cylinder  is  produced  by  pressurizing  the 
driver  section  of  the  shock  tube  shown  in  Figure  1  and  then  explosively 
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perforating  the  steel  diaphragm  at  the  interface  of  v..ie  dr  iver  section 
and  the  expansion  section  (Item  A  in  Figure  1).  Despite  the  fact  that  the 
thermal  source  was  located  downstream  from  the  target,  some  degradation  of 
the  leading  spike  of  the  shock  wave  occurred  due  to  the  presence  of  heated 
air  in  the  vicinity  of  the  cylinder. 

The  following  measurements  were  made  on  the  cylinder  during  the  tests 
pressure,  axial  and  circumferential  strain,  temperature,  and  thermal  flux. 
In  addition,  two  high  speed  cameras  were  used  to  monitor  the  sequence  of 
events  occurring  ir  the  simulator.  Pre-  and  post-shot  measurements  of  the 
radius  of  the  cylinders  were  made  at  384  points  on  each  cylinder.  More 
detailed  descriptions  of  the  instrumentation  are  contained  in  Reference  3. 

The  cylindrical  targets  were  made  from  1.016  mm  thick  sheets  of 
aluminum  6061-T6.  The  sheets  were  formed  into  cylinders  and  seam  welded 
on  the  side  farthest  from  the  loading  side.  The  inside  diameter  of  the 
cylinders  was  30.48  cm  (12  inches)  and  the  length  between  upper  and  lower 
clamps  was  80  cm.  In  order  to  make  the  ADINA  finite  element  modeling  of 
the  experiments  easier,  the  upper  and  lower  ends  of  the  cylinders  were 
clamped  securely  in  heavy  steel  end  discs.  The  end  discs  were  connected 
by  shafts  to  bearings.  The  bearings  were  mounted  in  two  heavy  bearing 
mounts.  To  minimize  the  disturbance  of  the  flow  past  the  cylinder,  the 
mounts  were  made  with  the  same  outer  diameter  as  the  target  cylinders. 

NUMERICAL  MODEL 

The  numerical  modeling  of  these  tests  was  done  with  the  ADINA  finite 
element  program  [4].  The  element  type  used  was  the  3-D  brick  element, 
which  has  a  thermo-elastic  formulation.  Symmetry  was  assumed  about  a 
plane  passing  through  the  cylinder  axis  oriented  in  the  direction  of  air- 
heat  flow.  Another  plane  of  symmetry  was  assumed  passing  orthogonal  to 
both  the  previous  plane  and  passing  through  the  Z-coordinate  L/2,  where  L 
is  the  active  length  of  the  cylinder.  With  these  assumptions,  only  one 
quarter  of  the  cylinder  was  required  for  the  finite  element  model.  Actu¬ 
ally,  there  was  approximately  a  10  per  cent  variation  in  heating  along  the 
axis  of  the  cylinder  due  to  nearby  reflecting  surfaces.  However,  this 
variation  was  Ignored  in  the  numerical  simulation. 

The  quarter  cylinder  was  modeled  with  72  brick  elements,  8  in  the 
axial  direction,  9  in  the  circumferential  direction,  and  1  in  the  thick¬ 
ness  direction.  The  elements  had  a  midside  node  on  each  edge  and  the 
number  of  integration  points  in  the  element  was  3  in  each  direction,  for 
a  total  of  27  per  element.  The  outside  surface  of  the  undeformed  mesh  is 
indicated  In  Figure  2  by  the  crosses  and  dotted  lines.  The  crosses 
represent  the  nodal  points  on  the  outside  surface  and  the  dotted  lines  are 
the  undeformed  mesh  boundaries. 
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Figure  2.  Mode  Shape  for  the  7th  Eigenfrequency  (3961  cycles/seconds) 

The  lower  eigenfrequencies  and  mode  shapes  were  evaluated  for  use  in 
a  mode  superposition  analysis  of  the  elastic  response  of  the  cold  cylinder. 
Knowing  the  periods  of  the  lower  eigenfrequencies  is  also  useful  in 
determining  the  approximate  time  required  for  a  peak  deflection  to  occur. 
The  periods  detex nined  ranged  from  3.66  ms  for  the  lowest  mode  to  1.47  ms 
for  the  eleventh  mode.  Shown  in  Figure  2  is  mode  number  seven,  which  has 
the  simpliest  geometric  shape  (egg-shaped  cross  section) . 

The  finite  element  calculations  for  the  dynamic  response  of  the 
cylinder  were  made  using  the  explicit  central  difference  time  marching 
scheme.  The  time  step  was  obtained  from  the  Courant  stability  condition. 


where  AJt,  is  the  distance  between  the  two  closest  nodes,  E^  is  Young' r 

modulus,  and  p  is  the  material  density.  For  the  higher  order  elements 
used  here,  we  found  that  this  time  increment  was  unstable  and  we  used  ?ne 
half  of  the  Courant  At,  which  was  stable.  This  time  increment  was  50  ns. 

The  constitutive  properties  of  the  cold  aluminum  6061-T6  material 
were  measured  in  an  Instron  testing  machine.  Tensile  test  specimens  were 
cut  from  the  same  sheet  of  material  used  for  fabrication  of  the  cylinders. 
Properties  were  measured  in  the  rolling  direction  and  perpendicular  to  the 
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rolling  direction.  The  average  of  these  tensile  properties  was  used  in  the 
ADINA  calculations.  Properties  measured  in  this  way  were  the  Young's  modu¬ 
lus,  bilinear  corner  yield  stress,  plastic  modulus  and  Poisson's  ratio. 


Dynamic  temperature  dependent  constitutive  properties  were  measured 
in  the  DNA  Tri-Service  Flash  Test  Facility.  In  these  tests,  samples  of  the 
material  were  heated  in  a  fraction  of  a  second  and  strained  through  the 
yield  point  several  seconds  later.  These  times  were  close  to  the  tempera¬ 
ture-time  sequence  experienced  in  the  thermal/blast  simulation  tests.  The 
yield  points  determined  in  these  tests  were  consistent  with  other  data 
[5],  However,  the  plastic  tangent  moduli  appeared  to  be  higher  than 
expected.  Due  to  this  result  and  to  difficulty  in  keeping  strain  gages 
bonded  at  the  high  temperatures,  we  chose  to  await  further  confirmation 
prior  to  using  the  data.  The  following  temperature  dependent  stress-strain 
data  have  been  derived  from  Reference  5.  No  temperature  dependent  data  for 


Temperature  Density 
°C  g/cm3 

28  2.703 
100  2.689 
150  2.678 
200  2.667 


Yield 

Young's  Modulus 

MPa 

GPa 

301.0 

64.73 

275.4 

62.20 

259.1 

60.80 

240.5 

59.29 

the  plastic  tangent  modulus  and  Poisson's  ratio  were  available.  The  cold 
properties  measured  in  the  Instron  testing  machine  of  650  MPa  and  0.3285, 
respectively,  were  used  for  all  temperature  ranges. 

The  ADINA  calculations  were  performed  using  a  bilinear  von  Mises 
kinematic  hardening  material  model.  This  model  required  the  specification 
of  the  materials  properties:  Young's  modulus,  Poisson's  ratio,  initial 
yield  stress,  and  plastic  tangent  modulus. 

RESULTS 

The  three  types  of  tests  and  analyses  will  be  discussed  in  the  same 
order  as  the  experiments  were  performed.  The  input  head-on  pressure  for 
the  blast-only  test  is  shown  in  Figure  3.  The  maximum  stagnation  over¬ 
pressure  was  60  kPa.  Shown  in  Figure  4  is  the  short  term  reflected  pres¬ 
sure  on  the  cylinder  at  zero  degrees.  The  peak  load  is  104  kFa  (15  psi) . 
Figure  5  shows  the  pressure  loading  on  the  side  of  the  cylinder  (90°). 
Loading  data  were  recorded  at  15°  Intervals  around  the  cylinder.  Since 
the  finite  element  calculations  required  loading  at  20°  intervals  around 
the  circumference  starting  at  10°,  these  data  were  interpolated  in  both 
circumferential  angle  and  time  to  supply  the  pressure  loading  for  ADINA. 
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Experimental  measurements  of  the  final  radial  deflection  of  the  cyl¬ 
inder  to  the  blast  load  were  made  and  the  results  are  shown  in  the  top  and 
middle  part  of  Figure  6  for  the  right  and  left  halves  of  the  cylinder.  The 
deflections  shown  have  been  magnified  by  a  factor  of  ten  to  make  the  de¬ 
formed  shape  more  apparent.  Deformation  of  the  actual  cylinder  was  barely 
noticeable  to  the  eye.  The  resulting  deformation  pattern  was  only 
approximately  symmetric  about  the  zero  degree  line  and  about  a  plane  at 
Z  »  L/2.  For  Figure  6,  the  loading  direction  is  from  the  top  for  all  three 
parts  of  the  figure.  The  deformed  shapes  for  the  three  subplots  at  the 
Z  »  L/2  center  cross  section  are  shown  on  the  left  side  of  Figure  6.  The 
ADINA  calculations  have  been  carried  out  to  2.0  ms  and  the  deflection 
pattern  at  that  time  is  shown  on  the  right  side  of  Figure  6.  Since  no 
plasticity  was  observed  in  this  calculation,  this  is  not  the  final  predict¬ 
ed  shape  of  the  cylinder,  but  is  approximately  the  maximum  deflection 
occurring  during  the  excitation  of  the  simulated  cylinder. 

In  the  second  test,  the  cylinder  was  exposed  to  a  thermal  load  only. 
Shown  in  Figure  7  is  the  flux  and  the  fluence  received  by  the  flux  gage 
positioned  off  to  the  side  of  the  cylinder.  The  solid  line  is  the  flux  as 
recorded  by  the  gage.  The  fluence  was  obtained  by  integrating  the  flux 
curve.  The  nominal  25  cal/cm2  fluence  on  the  target  was  obtained  by  re¬ 
lating  the  fluence  to  previous  calibration  tests.  Figures  8  and  9  show  the 
temperature  rise  indicated  by  thermocouples  at  Z  »  L/2,  0  ■  0"  and  90°  re¬ 
spectively.  The  nominal  temperature  rise  was  225°C  at  0°  and  A5°C  at  90°. 

The  ADINA  model  for  the  thermal-only  loading  was  the  same  mesh  config¬ 
uration  as  for  the  blast-only  loading.  However,  instead  of  the  pressure 
load  applied  to  the  exterior  of  the  cylinder,  an  initial  displacement  of 
the  Z-coordinates  of  the  cylinder  proportional  to  the  temperature  increase 
was  assumed.  A  prescribed  displacement  was  applied  to  the  symmetry  end  of 
the  cylinder  at  Z  *  L/2.  This  artifice  was  used  in  lieu  of  performing  an 
ADINAT  [6]  calculation  to  determine  the  temperature  in  the  cylinder  as  a 
function  of  time.  Earlier  calculations  using  ADINAT  showed  that  the  heat 
penetrates  through  the  thickness  of  1.016  mm  in  a  time  which  is  short 
compared  to  the  structural  response  time.  Based  on  this  result,  tempera¬ 
ture  was  assumed  to  reach  equilibrium  through  the  thickness  instanta¬ 
neously.  Initial  displacements  were  calculated  for  each  nodal  point  based 
on  the  fluence  received  at  that  point.  The  displacements  were  calculated 
relative  to  the  fixed  end  according  to  the  formula 

00 

w(e,z)  -  —  J  *(e,t)  dt, 

o 

where  a  is  the  coefficient  of  thermal  expansion,  h  is  the  thickness,  p  is 
the  density,  C  is  the  heat  capacity,  and  $(0,t)  is  the  incident  flux. 
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Measurements  on  the  post-test  thermally  loaded  cylinder  showed  a 
slight  degree  of  buckling  near  the  ends  of  the  cylinder  on  the  near  side. 
The  deformation  pattern  is  shown  in  the  top  two  subplots  of  Figure  10  and 
is  nearly  symmetric  in  both  the  Z  and  0  coordinates.  The  displacements 
from  the  original  geometric  positions  are  magnified  by  10X  in  Figure  10. 

The  deflection  was  small  and  was  most  apparent  to  the  unaided  eye  near  each 
end  where  the  largest  changes  in  slope  occurred.  The  deformation  cross 
sections  at  the  mid-point,  Z  *  L/2,  are  shown  on  the  left  side  of  Figure  10. 
The  displacement  was  outward  at  0  -  0°  and  inward  at  0  -  90°. 

The  ADINA  simulation  of  the  thermally  loaded  cylinder  was  carried  out 
to  1,1  ms.  Plasticity  occurred  early  in  the  calculation  and  peak  deflec¬ 
tions  were  observed  at  approximately  0.725  ms.  The  deformed  shape  at  this 
time  as  predicted  by  ADINA  is  shown  in  the  lower  part  of  Figure  10.  The 
apparent  cut  shown  in  Figure  10  at  Z  *  L/2  is  the  prescribed  constant 
Z-displacement  which  simulates  the  thermal  expansion.  The  prescribed  dis¬ 
placement  ranged  from  1.8  nun  at  0  »  0*  t  ,  zero  at  0  =  120°.  The  deforma¬ 
tion  cross  section  at  Z  *  L/2  showed  a  similar  trend  as  did  the  experi¬ 
mental  results.  However,  no  buckling  was  predicted  by  ADINA  near  the 
Z  -  L/20  point,  contrary  to  what  was  observed  in  the  experiment.  An  incip¬ 
ient  buckling  deflection  of  approximately  twice  the  cylinder  thickness 
appeared  to  be  forming  at  Z  ■  L/2,  0  «  100°  in  the  ADINA  results. 

In  the  final  test  in  the  BRL  Thermal/Blast  Simulator,  the  cylinder 
was  exposed  to  the  thermal  radiation,  rotated  and  exposed  to  the  air  shock. 
The  reference  flux  and  temperature  histories  at  0°  and  90°  are  shown  in 
Figures  11-13.  The  long  term  stagnation  pressure  record  is  shown  in 
Figure  14.  The  pressure  loadings  on  the  cylinder  at  0  -  0“  and  90°  are 
shown  in  Figure  15  and  16.  Comparing  Figures  15  and  16  with  Figures  4  and 
5,  one  notices  that  the  pressure  load  rises  much  more  gradually  in  the 
presence  of  the  heated  air  surrounding  the  cylinder.  From  Figures  14  and 
3,  one  sees  that  the  long  term  pressure  produced  by  the  shock  tube  is  only 
slightly  less  for  the  thermal/blast  test  than  was  the  case  for  the  blast- 
only  test.  As  shown  in  Figures  12  and  8,  the  peak  temperature  rise  in 
thermal/blast  test  was  170aC  and  in  the  thermal-only  test  it  was  225°C- 

The  final  deformed  shape  of  the  cylinder  after  the  combined  thermal/ 
blast  test  is  shown  in  Figure  17.  The  deformation  is  symmetric  with  re¬ 
spect  to  the  two  symmetry  planes.  The  deflection  magnification  in  Figure 
17  is  only  5X,  one  half  of  the  magnification  used  in  Figures  6  and  10.  At 
Z  -  L/2,  the  peak  radial  deflection  for  the  thermal/blr.st  test  is  approxi¬ 
mately  4  times  as  large  as  for  the  blast-only  deflection  and  10  times  as 
large  as  the  thermal-only  deflection.  ADINA  calculations  for  the  thermal/ 
blast  simulation  have  not  been  made  and  no  comparisons  are  available. 


;  0 


GREGORY  &  PEARSON 


Figure  12.  Thermal/Blast  Test,  Temperature  History  at  Z  «  L/10,  0 

TEST  THB9WL/BUST 


TDC.S 


Figure  13.  Thermal/Blast  Test,  Temperate  e  History  at  Z  ■  L/2,  6 
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TEST  THEBHAL/SUTT 


Figure  14.  Thermal/ Blast  Test,  Head-On  Pressure 

TEST.  TKEJWAL/BUST 


Figure  15.  Thermal/Blast  Test,  Diffraction  Phase  Loading  at  0 


TEST:  THERMAL  JUST 


Figure  16.  Thermal/Blast  Test,  Diffraction  Phase  Loading  at  90 
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LOADING  DIRECTION 

i 


PRE-TEST  TARGET  SHAPE  -  RIGHT  HALF 


POST-TEST  TARGET  SHAPE  -  RIGHT  HALF 


POST-TEST  TARGET  SHAPE  -  LEFT  HALF 

DEFLECTION  MAOSIIFICATION  SX 

Figure  17.  Comparison  or  Pre-Test  and  Post-Test  Shapes  for 
Thermal/Blast  Test 
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CONCLUSIONS 

The  series  of  tests  showed  that  the  effect  of  combined  thermal  and 
blast  loadings  at  realistic  time  separations  produces  a  nonlinear  response 
in  thermally  thin  targets.  Moreover,  the  shape  of  the  deformed  cylinders 
is  markedly  different  between  the  three  tests.  The  heated  cylinder  tended 
to  displace  outward  on  the  near  and  far  sides;  whereas,  the  combined  test 
resulted  in  a  dished-in  surface  on  the  near  side. 

The  results  obtained  thus  far  with  ADINA  are  not  as  satisfying  as  was 
originally  expected.  However,  the  permanent  deflections  produced  in  both 
the  tnermal-only  and  blast-only  tests  were  very  small.  Furthermore,  the 
general  shape  of  the  deformations  predicted  by  ADINA  was  correct.  It  is 
evident  that  we  need  to  develop  a  Kirchhoff  shell  element  whose  stable 
time  step  is  not  limited  by  the  rotatory  inertia.  This  would  improve  the 
efficiency  of  the  ADINA  calculation  significantly. 
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